Large radius exciton in single-walled carbon nanotubes 



Vadym M. Adamyan* Oleksii A. Smyrnovt 
Department of Theoretical Physics, Odessa I. I. Mechnikov National University, 
2 Dvoryanskaya St., Odessa 65026, Ukraine 

June 30, 2007 



Abstract 

The spectrum of large radius exciton in an individual semiconducting single-walled carbon nanotube 
(SWCNT) is described within the framework of elementary potential model, in which exciton is modeled 
as bound state of two oppositely charged quasi-particles confined on the tube surface. Due to the parity 
of the interaction potential the exciton states split into the odd and even series. It is shown that for the 
bare and screened Coulomb electron-hole (e-h) potentials the binding energy of even excitons in the ground 
state well exceeds the energy gap. The factors preventing the collapse of single-electron states in isolated 
semiconducting SWCNTs are discussed. 

PACS number(s): 78.67. Ch 

1 Introduction 

Many experimental papers on the optical absorption in SWCNTs describe obtained results in terms of band-to- 
band direct transitions between single particle states, though it is clear that the inherent to ID systems strong 
interparticle interaction cannot be neglected. It seems obvious that the strong electron-hole attraction should 
bind electron-hole pairs in SWCNTs into Wannier-Mott like excitons. Moreover, the exciton contributions were 
already revealed experimentally in optical absorption spectra pQ, [2], and in spectra of fluorescence [3]- [5] of 
individual SWCNTs, as well the exciton properties were studied by the Raman spectroscopy [BJ. There are also 
some works devoted to the theoretical study of excitons in CNTs [7]-[T2]. However, as it follows from results of 
the latter a simple translation of basic hydrogen-like models of 3D large radius excitons fails to function in one 
dimension without a certain specification. Remind that once the centrum of mass has been removed and the 
screening effect from the tube charges is ignored the ID model exciton Hamiltonian may be formally given by 
the expression 

H - ------ 

2 M dz 2 e \zY (1 ' 

where /j, is the exciton reduced effective mass and e is the dielectric constant of medium surrounding the tube. 
The functions <)>{z) from the domain of differential expression (jl.ip in l/2(— oo, oo) are twice differentiable on 
the semi-axis (— oo,0), (0, oo), and belong together with their derivatives for each 5 > to the subspaces 
L/2(— oo,— 8), L2(5,oo), and also satisfy a certain boundary conditions at the point z = 0, which provide the 
self-adjointncss of H . If there are no reasons for breaking of the tube reflection symmetry or, in other words, 
of the action and reaction law, then only those self-adjoint extensions of the above differential operator are 
physically admissible, for which the subspaces of even and odd functions from Li2(— 00,00) are invariant with 
respect to H. 

For less singular than the bare Coulomb even potentials V(z) = V{\z\), for example, for potentials satisfying 
the condition 



J \V{z)\ 2 dz < 00, < L < 00, 
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the corresponding Hamiltonians H could be represented as the sum of the standardly defined self-adjoint 
operator of kinetic energy and the subordinated operator of potential energy that is the multiplication operator 
by V(z). In such cases the functions from the domains T>(H) of H are continuous and have continuous first 
derivative at z = 0, and the corresponding odd and even functions from T>(H) should satisfy the natural 
boundary conditions <^>(0) = and <fi'(0) = 0, respectively. The spectra of bound states are then the set of 
negative eigenvalues of the boundary problem for the differential equation 

on the semi-axis (0, oo) with the boundary conditions cf>(0) = for the odd series and (f>'(0) — for the even 
series, respectively. 

For the Coulomb potential in one dimension the operator of potential energy is not subordinated to that of 
kinetic energy and the Hamiltonian H cannot be simply represented as their sum, and hence the Hamiltonian 
in this case becomes indeterminate. Among self-adjoint extensions of the differential operator (jl.ip there are 
infinitely many of those for which subspaces of odd and even functions are invariant. These extensions differ 
in self-adjoint boundary conditions at z = and, accordingly, they can be distinguished by energies of their 
ground states. Since functions from the domain of any such extension are non-differentiable at z = 0, it is 
impossible without additional physical considerations to single out the unique "correct" among suitable self- 
adjoint boundary conditions at z = 0. 

For the bare Coulomb potential one of possible extensions compatible with z-inversion symmetry is the 
decaying extension Hq defined by the boundary condition: (f)(0) = 0, that is for any <f> from the domain of Hq 
we have 

(flo^) (*) = -(h 2 /2fi)4>"(z) + (e 2 /ez)(f,(z), z < 0; 

(Sot) (z) = -(h 2 /2vW(z) - (e 2 /ez)(f>(z), z > 0; ( L2 ) 

<K-o) = <K+o) = o. 

This Hamiltonian was obtained in [13j . [14] as a result of formal passage to the limit for some sequences 
of ID Hamiltonians with regularized at z — Coulomb potentials. Note, that not only subspaces of odd and 
even functions, but also the subspaces of functions with supports on the positive and negative semi-axis are 
invariant with respect to Hq. In other words, for the zero boundary condition at the origin Hq is isomorphic 
to the orthogonal sum of two reduced "classic" Schrodinger operators for s-states of the hydrogen atom. As 
follows, the negative spectrum for this Hamiltonian is the Balmer series, each eigenvalue of which is doubly 
degenerate. It is worth mentioning, that if the states of electron-hole pair would be governed by the Hamiltonian 
Hq, then excitons in a tube would be subdivided into conserved "left" and "right" ones subject to the positional 
relationship of electron versus hole. 

In Section 2 of this paper we consider another extension Hi , which in the odd sector coincides with Hq but 
in the even sector is defined on the subset of continuous functions satisfying at z = the boundary condition 

lim — [(1 - 2Az\n{2A\z\)) Mz)] 

A d ^ (1-3) 

= lim— \(l + 2Azln(2Az))cj)(z)} = 0; A = e 2 [i/h 2 . 

z[Q dz 

For Hi the spectrum of bound states of the even series appeared to be close to that for the two-dimensional 
hydrogen atom [T5j for the states with zero angular momentum. This fact as well as the transition of (| 1 . 3[) 
to the Neumann condition (/>'(— 0) = <fi'(+0) = as A — > is not yet a valid reason to consider Hi as an 
appropriate primordial Hamiltonian for the large radius exciton in nanotubes. However, the modified electron- 
hole interaction potential V(z), that accounts, that these particles actually are not pointwisc and their charges 
are smeared along infinitesimal narrow bands on the tube surface, appeared to be locally quadratically integrable. 
In the case of nanotubes of small diameters the Hamiltonian with this potential gives the energies of ground 
and first excited states of the standardly defined even series, which differ slightly from those for H\. 

However, it turned out that the ground state energy of even excitons, calculated for individual semiconducting 
carbon nanotubes in vacuum with this potential and without account of the effect of screening by the nanotube 
electrons, are just two times greater of the energy gaps. Therefore, in Sections 3 and 4 we consider different 
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forms of screening of the electron- hole interaction inside individual semiconducting nanotubes (e.g., we calculate 
the dielectric function of individual SWCNTs). Results on the ground state of even excitons, given in Section 5 
for some individual carbon nanotubes, show that the account of screening does not help and the binding energy 
of even excitons remains greater of the energy gap. This may mean the instability of single-electron states 
in isolated semiconducting carbon nanotubes in the vicinity of the energy gap against the exciton formation. 
In the last section of the paper we discuss factors preventing the collapse of single-electron states in isolated 
semiconducting SWCNTs. 



2 Exciton spectrum and eigenf unctions in the Coulomb limit 

Let ijj v (k,r) and ip c (k,r) be the Bloch wave functions of the valence and conduction band electrons of a 
semiconducting nanotube, respectively. Remind that 

ip Vt c(k, r) = exp(i/cz)u v , c (fc, r), 

where u vc (k, r) are periodic functions with the period a along the tube axis, which is assumed to coincide with 
the z-axis. The wave functions of rest exciton can be represented as the following superposition: 

7r j a 

*(n,r 2 )= J $(fc)V> c *(fc,n)Vv(fc,r 2 )dfc. (2.1) 

— 7T / a 

The envelope function 3>(fc) in (|2.ip satisfies the equation 

7T / a 

{e c (k) - e v (k))$(k) + I J(k,k')$(k')dk' = E eKC $(k), (2.2) 

— tt / a 

where e v (k) and e c (k) are band energies of electrons with quasi- momentum k and the kernel 

J(k,k')=- Urn y / / Vc(fc,ri)VC(fc,r 2 ) , - r < (k\ n)^, r 2 )dndr 2 , 
l^oo L J J \n - r 2 | 

Eg = E 2 x (0 < z < L), 

corresponds to the two-particle system interacting with each other through the bare Coulomb potential. 
In the so-called long-wave approximation (|2.2p takes the form 

2 °° 

(^ S + $(fc) + h j V(k-k'Wk')dk' = E cx Mk), (2.3) 

— oo 

where E g and /i are the gap width and the reduced effective mass of electron and hole, respectively, and V(k) 
is the Fourier transform of the effective potential 

V(z) = 



1/2 

((xi ~ x 2 ) 2 + (yi - y 2 ) 2 + {z + zi- z 2 f) i 2 4) 

x |u c (0 ) ri)| a |u v (0,r 2 )| 2 dridra. 
We see from p.3p that ()2.2p is equivalent to ID Schrodinger equation 

h 2 

- ^"i z ) + V{z)^{z) = £cf>(z), £ = E cxc - E B , (2.5) 

on the real axis. Note, that independently on the tube radius and chirality 

V(z)\ ,^- e y\z\ + o(l/\z\). 
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Contrary to 3D case [H] direct using of the equation (|2.5|) with potential Vq(z) = —e 2 /\z\ for modelling of 
exciton states in nanotubes is impossible without a more accurate definition of the exciton Hamiltonian for 
short distances between electron and hole. The matter is that due to the Coulomb singularity of Vq(z) the 
one-dimensional Schrodinger operator for a hydrogen-like system remains indeterminate without imposing of 
certain (self-adjoint) boundary conditions onto wave functions at the point z = 0. So to define the ID exciton 
Hamiltonian we should either specify such a boundary condition or " soften" the singularity of potential at short 
distances with respect to the expression (|2.4j) and screening effects from the tube's electrons. As it was mentioned 
above, for parity of the Coulomb potential the exciton states are split into two series: even </>(— z) = <p(z) and 
odd 4>(—z) — —<f>(z). Despite of the Coulomb singularity at z = 0, any solution of the equation 

+ -r-r - K 2 U = 0, 
dz V M / (2.6) 

k = fie 2 1 nh 2 , k — \J 2/i\£\/h 2 

has continuous left and right limits at z = 0. Therefore continuous solutions of the odd series must satisfy the 
boundary condition: 

0(0) = 0. (2.7) 

Thus, the spectrum of bound states of the odd series coincides with that of the bound s-states of a hydrogen-like 
atom and this is also true for the corresponding wave functions on the positive semi-axis (up to the factor 1 / \/2) ■ 

However, with the Coulomb singularity the part of Hamiltonian on the subspace of even functions is not 
uniquely determined even under the condition, that the functions from its domain are continuous everywhere 
on the real axis, including the point z = 0. If the potential in the concerned problem would be nonsingular or, 
at least, integrable, then the boundary condition <fi'(0) — would be natural for the determination of this part. 
However, any non-zero at z — solution of (|2.6[) is not differentiable at z = 0. Q 

Attempts to choose a proper boundary condition at z — for the even part of Hamiltonian by some 
parity-preserving regularization of the Coulomb potential in the S- vicinity of the origin give non-unique results, 
depending on ways of regularization and passing to the limit as 5 — > 0. To see this let us consider the regularized 
potential 

' -e 2 /\z\, \z\>6; 
-e 2 /\S\, \z\<S. 



V S (z) 



Since Vq(z) < Vs(z), then the least eigenvalue £ Q {8) of the odd series for the Schrodinger operator Hs with the 
potential is not less than the energy of the ground state of the hydrogen-like atom, that is 

Let us define further the even part of the Schrodinger operator Hs with potential V$, assuming that even 
functions from its domain satisfy the boundary condition: 

0'(-O) + M>(-o) = -0'(+o) + M(+o) = 0. (2.8) 

Taking any 

we can arrange by a suitable choice of hs in (|2.8p that £ c be an eigenvalue of Hs- To this end we note, that for 
z > 5 the eigenfunction 4> c (z), corresponding to the eigenvalue £ a , coincides up to a constant factor with the 
decreasing as z — > oo solution of (|2.6[) with £ replaced by £ c , that is 

cf> e (z) = C W keA/2 {2K e z), 

where is the Whittaker function. At the same time for < z < S we have 

<p e (z) = C [cos qz + (h s /q) sin qz] , 



q=sj{2n/h 2 ) (e 2 /6-\£ e \). 



1 Actually, in one dimension the representation of Hamiltonian as the sum of operators of kinetic and potential energies is strictly 
speaking impossible for potentials with Coulomb singularities. 
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The continuity condition for the logarithmic derivative of (freiz) at z = S yields 



h s = q 



qW kBtl / 2 (2K e 6) cosqS - 2K e W^ 1/2 {2K e 8) sing6 
2n e k e ln(2«; c fc c (5). 



s | o 



(2.9) 



As £ c < £ (5) and eigenvalues of even and odd series alternate we conclude that £ e is the least eigenvalue of Hg. 
We see, that by an appropriate choice of hs we obtain a sequence of Hamiltonians with regularized potentials 
and the fixed least eigenvalue. 

As a nearest analogue of the boundary condition </>'(0) = for wave functions of the even series we take with 
account of (|2.9p : 

lim — [(1 + 2Az ]n(2Az)) Mz)] = 0, (2.10) 
z^o dz 

where A = e 2 fi/h 2 . The Schrodinger differential operator H in L 2 (0,oo) defined by this boundary condition is 
self- adjoint. 

Indeed, using the von Neumann formulae [17j it is easy to verify, that the functions from the domain T>(H) 
of each self-adjoint extension H of the differential operator (jl.ip in 1/2(0, oo) are continuous on the semi-axis 
[0,oo). 

Let 4>i,4>2 € V{H), that is H4>i,H<j> 2 € L 2 (0, oo). By continuity of <j>\,<j)2 at z = and (|2.10[) we get 



(H Vl ) (z) V * 2 (z) - <pi(z) (H V2 )* 



dz 



lim 

510 



(H Vl ) (z)p* 2 (z) - tptiz) (H^y (z) 



dz 



= lim [(pi(z)^ 2 *(z) - ip'^LpHz)]^ 



= lim 



1 



sio 11 + 2Az \n2 Az 



<Pi(z)l-[(l + 2Az\n2Az)<p%(z)] 
dz ' 



-<P2(z) — [(l + 2Azln2Az)ipi(z)} 
dz 



0. 



Therefore if is a symmetric operator. Let us assume that H is not self-adjoint. Then for each non-real u 
there is a solution of the equation 

h 2 d 2 W e 2 

W = ujW. 

2n dz 2 \z\ 

which belongs to 1*2(0, oo) and orthogonal to the linear set — iSj T>{H) [17]. But it is easy to verify as above, 

that such solution is identically equal to zero. Thus energy levels of the even series are defined as eigenvalues 
of Hamiltonian 

H-- — —~ — 
2fidz 2 \z\' 

on set of twice differentiable functions <j>(z) at semi-axis (0,oo), which satisfy boundary condition (|2.10|) . At 
semi-axis (— oo, 0) we, naturally, consider even continuation of corresponding eigcnfunctions. 

Evidently, the above choice of the even part of the Hamiltonian is not exceptional. As it was mentioned 
above, for example, we can take the even extension onto the negative semi-axis of the wave functions satisfying 
the zero boundary condition at z = and obtain in this way an even part of Hamiltonian with the same 
spectrum as that for the odd part [13] . The choice of a concrete boundary condition can be done exceptionally 
on the base of physical reasons. 

Using the asymptotic expansion of the Whittaker function Wk y i/ 2 (^Kz) for z — > we get from condition (|2.7|) 
the eigenvalues of the odd series: 

— !— = 0; 1 - k = -n, n 0.1.2...: 
( } , , (2.11) 



Me 4 1 
2h n 2 
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and from condition (|2.10[) the eigenvalues of the even series: 



where p, according to (|2.10p , is defined from equation: 

00 1 1 

-£>y^ — r + 7; hip + 7-1 = 0, 

f^j(j-p) 2P (2.13) 

^p = n+ 1/2 + A(n), n = 0,1,2.... 

Here 7 ~ 0.5772 is Euler's constant, and A(n) - is slowly increasing function of the integer number n, which 
over the range n € [0, 10] obtains values from —0.013 to 0.0156. 

Corresponding normalized wave functions <j) n (z), which satisfy the equation (|2.6p and the condition (12. 7p 
are given by: 

1/2 / 2Az\ 2Az ( Az 



^ ( * )= U(n-l)»(n-l)!» J Ll ^\~)~ eX? \~) n = 2 ' 3 ' 4 "- (2.14) 
tf>i(z) = V2A 3/2 zcxp(-Az), 

where L] l _ 1 (2Az/n) is the generalized Laguerre polynomial. For the even series, according to (|2.6|) and ([2.10fi . 
we obtain: 

/ 2Az\ 

cf) p (z) = C p W Ptl/2 f — J , (2.15) 

where C p is normalization factor. 

The analytic simplification of (|2.4[) . which depends on the tube radius Rq but independent of its chirality is 
the potential 

vM = -j4n 11 — m> ^ 

that was obtained from (|2.4[) under the assumption that the charges of electron and hole participating in the 
formation of exciton are smeared uniformly along infinitesimal narrow bands on the tube wall. This potential is 
the simplest approximation to the bare Coulomb potential, which accounts the finiteness of the tube diameter. 
Note, that contrary to the bare Coulomb potential this one has only logarithmic singularity at the origin. Since 
Vro i z ) is an integrable function then solutions of the Schrodinger equation with this potential are continuously 
diffcrentiable at z = 0, and the boundary condition for the even series in this case is </>'(0) = 0. For nanotubes 
with rather small diameters the negative eigenvalues of equation (|2.5p with potential Vr (z) appeared to be close 
to those for equation (12. 6ft (see Section 5, table 1 and table 2). For the both of equations the minimal eigenvalue 
of the even series well exceeds the energy gap. This may mean that the single-electron states in semiconducting 
SWCNTs in the vicinity of the energy gap are unstable with regard to the formation of excitons. However, the 
screening of e-h interaction by the tube electrons could result in the shift of exciton levels into the gap. To 
make clear whether it is so we consider further different forms of screening of the potential (|2.16p . 



3 Nanotube dielectric function 

First we obtain the nanotube dielectric function within the framework of the Lindhard method (the so-called 
RPA), then in the limiting case of small wavenumber values we get the Thomas-Fermi screening theory for 
charged particles in semiconducting SWCNTs. 

Following the Lindhard method, to obtain the e-h interaction potential f{r), screened by the electrons 
of quasione-dimensional nanotube lattice, we consider the one-dimensional Fourier transform of the Poisson 
equation: 

(q 2 -A 2B )rtq,r 2D )=47r(p™ t (q,r 2D )+p iad (q,r 2D )), (3.1) 

where r 2 D is the transverse component of the radius- vector, q is the longitudinal component of wave vector, 
p ext {q, J"2d) is the one-dimensional Fourier transform along the tube axis of the density of extraneous charge 
p ext (z, t 2 d) and p lnd (q, J" 2 d) is that of the charge density induced by the extraneous charge. Further we will 
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assume that p ext is axial symmetric, p cxt (q, T2d) = p cxt (q, »"2d), and localized in the small vicinity of the tube 
wall. As follows ip(q, r 2 v) and p lnd (q, v 2 d) depend on r 2 v only through r 2 D and besides whatever the case p lnd 
is localized at the tube wall. By (|3.1|) the screened e-h interaction potential may be written as: 



r 2D ) = 4vr | (p cxt (g, r 2D ) + p ind (g, r 2D )) G (q, r 2D , r' 2D )dr' 2D , 



(3.2) 



where Go(q, r 2 u, ^ 2 d) = (l/27r)jKo(|?j 1t*2D — ^dI) ^ s the Green function of the 2D Hclmholtz equation, and Kq 
is the modified Bessel function of the second kind. 

Let E®(k) and ^ s (t*) = (1/^/N) exp(ifcz)u^ s (r) be the band energies and corresponding Bloch wave 
functions of the nanotube 7r-electrons and E s (k), 5'fe. s (r) be those in the presence of the extraneous charge. 
Then 



p ind (<Z,r 2D ) - -e I exp(-igz)^ [/(£ s (fc))|* fe , s (r)| 2 



(3.3) 



f{E»{k))m s {r)\*\ dz, 



where / is the Fermi-Dirac function, L is the length of CNT and s numbers single-electron bands (N is the 
number of unit cells in the nanotube). In the linear in (p approximation we get: 



(f_ B s ^(k,k-q, a) 



(3.4) 



x / M*. fc _ g s (z, r2u)u c -k,s'(z, r 2 u)dzip(q, R ), 
o 

where a is the longitudinal period of nanotube and 

u* k ,(z, r2T>)u v; k-q,s(z, r 2 D)dzdr 2 D, 



B s . s >(k, k- q,a) = 

e 2 o 

E z . tStS ,{k) = E c;s ,(k) - E v;s (k). 

Taking into account the axial symmetry of p ext (q, r 2 u) and p md (q, t*2d) and their localization near the nanotube 
wall (r2D = Rq) we obtain from (|3.2p and (|3.4[) that 



<p(q,R ) = $(q,R ) + 2I (\q\R )K (\q\R ) / p md (g, r 2 D)dr 2 D 



(3.5) 



where ip(q,Ro) is the Fourier transform of the electrostatic potential induced by p ext , and Iq is the modified 
Bessel function of the first kind. We see that: 



¥>(<7,-Ro) 



£-Ro,a(<7) ' 



7T / a 



— Tv/a 



\B s ^(k,k- q, a) 
Eg- SyS ' (k) 



-dk I (\q\R )K (\q\Ro). 



In the limiting case of small wavenumbers: 



\B s ^(k,k- q,a)\ 2 Q ps \U s>s >(k, a)\ 2 q 2 , 



|£/ s ,y(M)l = 



d 

u *c;k,s<{ z -, r 2 v)-7^;u v . k ^(z, r 2D )d2;dr2D 



e 2 o 



(3.6) 



(3.7) 
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Note that C/ SjS '(fc,a) is nonzero only for the mirror bands, that is U s , s >(k,a) = U s (k, a)S s ^ s i . Using the orthogo- 
nality of the Bloch wave functions and applying the Schrodinger equation for ^k,s(z, T2d) yields 



U.(k,a) 



ih 2 N 



m Q E g;StS (k) 



d 

** fe s (z, r 2 D)^*v ; fe, ;5 (^, t 2 d) dzdr 2 D- 



(3.8) 



E 2 



Hence, the screened quasione-dimensional electrostatic potential induced by a charge eo, distributed with the 
density: 



p e * x (r) = ^-6(z)5(r 2 n-R ), 
in accordance with (|3.6[) and (|3 . T[) . is given by the expression 



<p(z) = 



e o f Io(\q\)K (\q\)exp(iqz/Ro) 



7tR J l+g a q 2 I (\q\)K (\q\) 



dq, 



(3.9) 



with 



e 2 h A x , 
9a = 9 n2 



7r j a 



nrri^R 2 



-7r / 'a 



d_ 

dz 



dft. 



(3.10) 



This potential was calculated using the single-electron energy spectrum and wave functions, obtained in [18j . 
The ground state exciton binding energy, calculated from (|2.5p with the screened potential (|3.9p . remains 
noticeably greater than the energy gap (see Section 5, table 4). 



4 Screening by free charges 

Free charges may appear in semiconducting nanotubes at rather high temperatures T. So here we will ob- 
tain the self-consistent screened potential of e-h interaction depending on the nanotube diameter and medium 
temperature. 

To take into account the screening of e-h interaction potential <p(r) by free charges (by intrinsic electrons 
and holes) we consider the Poisson equation: 

- A<p + K 2 R S(r - R )tp = 47re<5(r - r ), (4.1) 

where we suppose again, that the screening particles (electrons and holes) and the screened e-h pair itself are 
localized at the surface of cylinder (nanotube's wall) with radius Rq. Here k 2 = (47re 2 no/fcBT)(l/7ri?§), and 



sj 2^/^ST/27rft) exp(-£ g /2fc B T) 



7i () 

is the one-dimensional analogue of particle concentration in the intrinsic semiconductors. We assume, that 
CNTs can be treated as such semiconductors. 

Equation (|4.ip (without factor 47re) can be represented in the equivalent form: 



G(r, r„) = G (r, r„) - n 2 R a / G (r, r')S(r' - Ra)G(r', r )dr', (4.2) 



where Go(r, ro) = 1/4tt\i — ro| is the Green function of the Poisson equation without screening (k = 0). After 
averaging over axial and radial components of the radius-vector and several Fourier transforms, we obtain the 
following one-dimensional screened e-h interaction potential ip(z): 

00 

/ s _ 1 f ^ (fc)exp(ifcz) 
^ Z) ~V^J l + (27r)3/2( Kj R )2^ o(fc)/ 4 7re ^ 




Figure 1: The e-h interaction potentials versus the electron-hole distance z for the zig-zag nanotube (28,0): 
dashed line - the screened one calculated by (|4.5|) for T ~ 550°K; black circles - the screened potential (|3.9|) ; 
solid line - the Coulomb unscreened averaged potential 16[) : dot-dashed line - the bare Coulomb potential 
from (HU) 



where ifio(k) is the Fourier transform of the average unscreened potential (|2.16[) : 



2r, 



~ ru\ 4ne f 1 f exp(-ifcz) 



8n*J yfcj^ \~ z \( 1 + {4 Rt/P)sm 2 (a/2)) 1/2 ' (44) 
47re 



(2tt) 3 / ; 



;J (\k\Ro)K (\k\Ro), 



where Io(\k\Ro) and Ko(\k\Ro) are the same modified Bessel functions of the first and the second kind respec- 
tively. Hence, the e-h interaction potential screened by free charges for any semiconducting SWCNT is given 
by: 

OO 

^ [Z) -nR a J l + (KR )n (\k\)K (\k\) ™- [ ^ } 

— oo 

The screened potential (|4.5p can be used for the calculation of the large radius exciton binding energies in the 
ground and excited states for the large-diameter SWCNTs at high temperatures. The ground state binding 
energy, calculated from (|2.5j) with the screened potential (|4.5| , remains greater than the energy gap (see Section 
5, table 5). 

To compare different obtained potentials, we produce figure 1, that shows them plotted point by point for 
the semiconducting (28,0) nanotube in comparison with the bare Coulomb potential. Figure 1 also shows that 
the mentioned above screened potentials slightly differ from the bare Coulomb potential when the distance 
between electron and hole is large. This fact justifies the bare Coulomb large radius exciton model given in the 
beginning. 



5 Calculation results. Screening influence 

Electronic structure of nanotubes, electron and hole effective masses and energy gap magnitudes were obtained 
in [TS] within the framework of the zero-range potential method for the Bloch wave functions [19 . Using those 
values of effective masses and energy gaps, we have calculated the unscreened and screened e-h interaction 
potentials and corresponding exciton binding energies of the ground and excited states, which either explicitly 
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or implicitly depend on parameters of concrete semiconducting SWCNT (chirality, radius, reduced effective 
mass, band gap magnitude) and the temperature of medium (in Section 4). Here we present results of these 
calculations. 

Numerically calculated values of the exciton binding energies of the ground and excited states according 
to (|2.6p are given in table 1. 



Table 1: Exciton binding energies according to (12. 6|) . 



Chirality 


E g) eV 


^0;even; 


£i ; odd, eV 


W;evem 


^2;odd, eV 


(7,0) 


1.3416 


-2.8343 


-0.6722 


-0.294 


-0.168 


(6,5) 


1.1017 


-2.9253 


-0.6938 


-0.3034 


-0.1734 


(28,0) 


0.3674 


-0.9484 


-0.2249 


-0.0983 


-0.0562 



These results unambiguously show that the binding energies in the even ground state for any of the selected 
semiconducting SWCNTs are much greater than the corresponding energy gaps in the bare Coulomb limit (|2.6p . 

Further, the numerically calculated values of exciton binding energies at the ground and excited states 
according to the wave equation with the potential (|2.16p are given in table 2. It can be seen from table 2, 
that the discrepancies with the analogous results in table 1 are more considerable for nanotubes with larger 
diameters, because the wave equation with the potential (|2.16l) tends to (|2.6[) if i?o — > 0. 

Table 2: Exciton binding energies according to the Schrodinger equation with the potential (|2.16|) 



Chirality 


2i?o, nm 


M);even; 


£i;odd) eV 




£2;odd, eV 


(7,0) 


0.548 


-2.7894 


-0.5365 


-0.2955 


-0.1507 


(6,5) 


0.7468 


-2.2503 


-0.5071 


-0.2829 


-0.1488 


(28,0) 


2.192 


-0.7567 


-0.1668 


-0.0928 


-0.0486 



We can see also from table 2 that the ground state binding energies are larger than the corresponding energy 
gaps even if the finiteness of nanotubes is taken into account. 

The next table (table 3) shows that the exciton radii are comparable with the corresponding nanotubes 
diameters, thus they much greater than the nanotube lattice parameter 0.142 nm. Therefore the large radius 
exciton theory methods are appropriate for the treatment of the SWCNTs exciton problem. 



Table 3: Exciton radii r n ~ n/2A in units of 2Rq 



Chirality 


f"0;evcn 


rijodd 


1; even 


7"2;odd 


(7,0) 


0.4759 


0.9772 


1.4776 


1.9545 


(6,5) 


0.3383 


0.6948 


1.05 


1.3895 


(28,0) 


0.3556 


0.73 


1.1 


1.46 



As illustration we have calculated the binding energies for the (28,0) zig-zag nanotube with account of the 
nanotube dielectric function (table 4). 

The data from table 4 obviously show that the screening by nanotube band electrons is not enough for the 
ground state exciton binding energy to be less than the energy gap. 

The exciton binding energies at the ground and excited states for the semiconducting (28,0) SWCNT cal- 
culated using the potential (|4.5|) for T — 550°K are listed in table 5. 

Note, that the screened potential (|4.5|) may be used either for the semiconducting SWCNTs with narrow 
band gap (as zig-zag (3n, 0) SWCNTs) or for the large-diameter nanotubes (small gaps) or(and) at rather high 
temperatures, because only under these conditions the linear concentration no of free charged particles provides 
a perceptible screening. At T — 550°K the (28,0) nanotube has approximately one free charged particle per 
micrometer of its length, but even at these conditions the screening by free charges of the e-h interaction 
potential is much stronger than the screening by the all bound electrons of semiconducting SWCNTs (compare 
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Table 4: Exciton binding energies for nanotube (28,0) according to (|3.9[) and (|3.10p 



Chirality 


9a 


^0; even 5 


£l:odd, eV 


^l;even) 


^2:odd, eV 


(28,0) 


0.6 


-0.6869 


-0.1799 


-0.0952 


-0.0501 



Table 5: Exciton binding energies for the nanotube (28,0) according to (|4.5[) for T = 550°K 



Chirality 


kRq 


^0;even; 


^i ; odd, eV 


^1; even 3 


^2;odd, eV 


(28,0) 


0.38 


-0.5549 


-0.0642 


-0.0294 


-0.0133 



table 4 and table 5). Nevertheless, as it follows from the same table 5, even in this case the ground state exciton 
binding energy still exceeds the energy gap. 



6 Discussion 

In the all above examples the binding energy of the ground state of even excitons in isolated SWCNTs appeared 
to be much greater than the corresponding band gaps even with account of some screening effects by tubes 
7r-electrons. This should mean that the single-electron states in SWCNTs are unstable at least in the vicinity 
of the energy gap with respect to formation of excitons. Such conclusion might seem doubtful though we came 
to it applying similar arguments as in the case of 3D large radius excitons. There are three reasons due to 
which a partial destruction of band electrons states in semiconducting SWCNTs in reality is either absent or 
inconspicuous. 

First of all the account of dynamical screening, that is the frequency dependence of dielectric function, may 
return the all exciton levels into the band gap. This was shown in [20], where calculations of the exciton binding 
energy with the static dielectric function yielded also the exciton binding energy exceeding the energy gap. At 
the same time the self-consistent calculation with frequency dependent dielectric function gave according to [20 
a universal ratio of the exciton binding energy to the energy gap depending only on the resonance integral 70 but 
not on the nanotube radius (it equals 0.87 if 70 = 2.7 eV). By [20] the exciton binding energy cannot be larger 
than the energy gap because of the singularity of the frequency-dependent dielectric function e{uj) at uj = E g /h 
for the frequencies, corresponding to the direct transitions between the van Hove points of the tube single- 
electron spectral density. However, actually this argument is true only if the exciton binding energy obtained 
without account of dynamical screening gets into a small vicinity of the energy of allowed transition between 
such points. This is because the frequency dependent SWCNT dielectric function may only then become rather 
great. Otherwise as it follows from results of [21] the effect of dynamical screening is too small and the exciton 
state with the binding energy much greater than the energy gap transforms into a long-living resonance in the 
continuous spectrum of electron-hole pairs with opposite quasi-momenta. 

The second reason is the so-called environmental effect. In experimental works [1]-[B] (which used the meth- 
ods described in |3J) investigated individual nanotubes were not in vacuum but encased in sodium dodecyl 
sulfate (SDS) cylindrical micelles disposed in D2O. Because of these SDS micelles, which provided a pure 
hydrocarbon environment around individual nanotubes, the high permittivity solvent D2O did not reach nan- 
otubes. However, the environment of hydrophobic hydrocarbon "tales" (— C12H25) of the SDS molecules has 
the permittivity greater than unity. Following the figure 1A from [3] we considered a simple model of a SWCNT 
in a dielectric environment: a hollow, narrow, infinite cylinder with radius Ro in a medium with the dielectric 
constant e and found the potential (|2.16|) screened by the medium within the framework of mentioned model 
under the assumption about axially symmetrical charge localization at nanotube's (here - cylinder's) wall. The 
corresponding ID screened potential (p(z) is given by: 

00 

( , e f I (\k\)K a (\k\)exp(ikz/R ) 

^ ' ttRq J [eKt(\k\)Io(\k\) + Ii(\k\)K (\k\)]\k\ [ °' > 

—00 

where ij(|fc|) and fC,(|fc|) are the modified Bessel functions of the order j of the first and the second kind, 
respectively. We don't know the exact value of dielectric constant of the pure medium, which is formed from 
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the hydrocarbon "tales" of the SDS molecules. But for estimates we take the dielectric constants of the 
substances, which are also formed from similar hydrocarbon "tales", e.g.: petroleum (e ~ 2.1) or dodecane 
(e ~ 2) at 293°K (this temperature is very close to that used in [Sj-[5]), or polyethylene (e ~ 2.2 — 2.4). Using 
the potentials (|2.16|) and (|6.ip with e varying in the interval 2 — 2.4 we have got that the ground state exciton 
binding energy in the nanotube (8,0) (the energy gap equals 1.415 eV |18j ) is 3.06 eV in vacuum while with 
account of the environment it runs the interval 1.33— 1.06 eV and hence gets into the corresponding energy gap 
and becomes close to those in [10] (about 0.86 — 1 eV), even without account of static and dynamical dielectric 
screening of the potential (|2.16|) by nanotube electrons. Remind that results on the (8,0) nanotube in [10] 
are in good agreement with those obtained in [5] by interpolation of experimental data for another species of 
nanotubes. 

Further, taking the (7,5) nanotube we compare our results with the corresponding experimental data from [BJ, 
where individual SWCNTs were isolated in surfactant micelles of SDS in D 2 like in [3]. Our calculations for the 
(7,5) nanotube in vacuum yield 2.12 eV as the ground state exciton binding energy, while for the same tube in 
the SDS environment the binding energy calculated using the potential Q6.1|) gets into the interval 0.90 — 0.71 eV 
(the band gap for the (7,5) tube is 1.01 eV 18]) depending on e varying from 2 to 2.4. The obtained binding 
energy value is not far from that of [6] ~ 0.62 eV even without the account of static and dynamical dielectric 
screening of the potential (|2.16|) by the nanotube electrons. There is a comparison of experimental data on 
the exciton binding energies in the work [6] with the corresponding theoretical results of [11] . These results 
are well agreed. But again, in the work [llj the interparticle potential includes screening parameter denoted 
as k = 2. Besides, it is asserted in [TT] that the assumption of similar Coulomb parameters for SWCNTs and 
phcnyl-based 7r-conjugated polymers, used in this work, gives smaller exciton binding energies for SWCNTs. 
All the results listed in table 1 - table 5 of our work are related only to SWCNTs in vacuum. So let us turn 
to the experimental work [22] which deals with optical properties (photoluminescence) of SWCNTs suspended 
in air (near- unit dielectric constant). As it follows from [25], the relative discrepancies between the optical 
transition energies obtained in [22] and those obtained in [5] are not significant (about several percents). This 
result could be expected, since according to the usual self-consistent field approximations the interaction of a 
7r-electron with other electrons of a nanotube should be substantially compensated in the ground state by the 
interaction with the nearest ions. Evidently, the effect of this compensation is not sensitive to an environmental 
screening. However, for excited states such as excitons, where electrons and holes are at distances of the order 
of tube diameter, the environmental effect can be strong. 

Note thirdly that with the advent of N excitons in the tube the additional screening effect, stipulated by a 
rather great polarizability of excitons in the longitudinal electric field, appears. The elementary estimates show 
that the corresponding adding to the dielectric constant is 

where £^ is the binding energy of even exciton in the ground state and L is the length of a tube. We see that 
in the case of N ~ 10 per 100 nm of nanotube length Ae ~ 1 and therefore the lowest exciton binding energy 
occurs already inside the energy gap. This blocks further conversions of single-electron states into excitons. The 
shift of the forbidden band edges due to the transformation of some single-electron states into excitons results 
in some enhancement of the energy gap. As follows the optical transition energy En should be blueshifted as 
in [22]. A coarse estimate of this shift using the elementary relation 

TiViV 2 

gives AEn/ En ~ 10~ 2 . If the exciton gas in tubes is unstable with respect to transition into a one-dimensional 
electron-hole plasma, then for the account of screening effect produced by this plasma we can use the results 
of Section 4. For example, for the (8,0) tube even ten charges per 100 nm of its length (~ 0.1% of 7r-electrons 
number) reduce the ground state exciton binding energy to 0.12 eV and thus block spontaneous transitions to 
the exciton states. 

Thus we may conclude that the ground state of 7r-electrons in semiconducting SWCNTs in vacuum is formed 
by band electrons filling all the levels up to a certain level below the gap together with some amount of two- 
particle even excitations, which can form either a rare gas of excitons or electron-hole plasma. The additional 
screening effect induced by the exciton gas (or the one-dimensional e-h plasma) blocks further partial destruction 
of single-electron states. The environmental effect may return the even exciton binding energies into the energy 
gap and thus may remove two-particle excitations from the ground state of 7r-electrons in SWCNTs. 
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